Decoherence of an n-qubit quantum memory 
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We analyze decoherence of a quantum register in the absence of non-local operations i.e. of n 
non-interacting qubits coupled to an environment. The problem is solved in terms of a sum rule 
which implies linear scaling in the number of qubits. Each term involves a single qubit and its 
entanglement with the remaining ones. Two conditions are essential: first decoherence must be 
small and second the coupling of different qubits must be uncorrelated in the interaction picture. 
We apply the result to a random matrix model, and illustrate its reach considering a GHZ state 
coupled to a spin bath. 
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Physical devices capable of storing faithfully a quan- 
tum state, i.e. quantum memories, are crucial for any 
quantum information task. While different types of sys- 
tems have been considered, most of the effort is con- 
centrated in manipulating qubits as they are essential 
for most quantum information tasks For a general 
quantum memory (QM) it is necessary to record, store 
and retrieve an arbitrary state. For quantum commu- 
nication and quantum computation initialization in the 
base state, single qubit manipulations, a two qubit gate 
(e.g. control-not), and single qubit measurements are suf- 
ficient. For arbitrary quantum memories recording and 
retrieving the state is still an experimental challenge but 
initializing a qubit register and individual measurement 
is mostly mastered. Faithful realization of two qubit 
gates together with better isolation from the environ- 
ment are the remaining obstacles for achieving fully op- 
erational quantum technology. However the huge effort 
done in the community has born some fruit 0, E3] • 

A number of papers discuss how two qubit quantum 
memories are affected by coupling to some environment. 
In our previous work we study the case of two qubits 
interacting either with a spin bath [4( or a random ma- 
trix environment 0, Q . In both cases the loss of coher- 
ence and the loss of internal entanglement of the pair was 
studied. Many other papers focus on the loss of internal 
entanglement only, due to the interaction with a thermal 
bath (see Q and references therein) . A study of decoher- 
ence of a two-qubit system as measured by entanglement 
fidelity is presented in [8}. There the effect of fluctua- 
tions of the parameters of the Hamiltonian is analyzed. 
In U3] an n qubit QM coupled to a bath of harmonic 
oscillators is studied. 

In the present letter we give analytic expressions for 
the decoherence of a QM during the storage time. Specif- 
ically we discuss a QM composed of a set of individual 
qubits interacting with an environment. The expressions 



given are based on previous knowledge of the decoher- 
ence of a single qubit entangled with some spectator, 
with which it does not interact. Their validity is limited 
to small decoherence, i.e. large purity of the QM. Note 
that the latter is no significant restriction due to the high 
fidelity requirements for the application of quantum er- 
ror correction codes. We assume that the entire system 
is subject to unitary time-evolution, and that decoher- 
ence comes about by entanglement between the central 
system and the environment. Because of its simple ana- 
lytical structure decoherence shall be measured in terms 
of purity IlJ of the central system. Spurious interac- 
tions inside the central system are neglected. We shall 
rely heavily on our recent studies of decoherence of two 
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qubits 

A further and critical assumption is the independence 
of the coupling in the interaction picture. This is justi- 
fied if the couplings are independent in the Schrodinger 
picture or if we have rapidly decaying correlations due 
to mixing properties of the environment in the spirit of 
[ljj ]. Physically the first would be more likely if we talk 
about qubits realized in different systems or degrees of 
freedom, while the second seems plausible for many typ- 
ical environments. 

We generalize the concept of a spectator configura- 
tion 0, S Sj: The central system consists of two non- 
interacting parts, one interacting with the environment 
and the other not. This configuration is non-trivial if the 
two parts of the central system are entangled. In our 
case the first part will be a single qubit and the second 
part (the spectator) will consist of all other qubits of the 
QM. Our central result is then a decomposition of the 
decoherence of the full QM, as measured by purity, in a 
sum of terms describing the decoherence of each qubit in 
a spectator configuration. Apart from the above made 
assumptions, this result does not depend on any partic- 
ular property of the environment or the coupling. Thus 
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it can be applied to a variety of models. 

We test successfully the results in a random matrix 
model both for coupling and bath 0] . The general re- 
lation is obtained in linear response approximation and 
leads to explicit analytic results if the spectral correla- 
tions of the bath are known. Finally, we perform numer- 
ical simulations for two and four qubits, interacting with 
a kicked Ising spin chain as an environment, 

which can be chosen mixing or integrable. 

The central system (our QM) is composed of n qubits. 
Thus its Hilbert space is H qm = &)™ =1 Hi, where Hi 
are the Hilbert spaces of the qubits (dimHi = 2). The 
Hilbert space of the environment is denoted by Ti c . The 
Hamiltonian reads as 



H = H, 



qm 



H C + XV 



(1) 



Here, H qm = X^iLi where Hi acts on Hi, whereas H e 
describes the dynamics of the environment and 



(i) 



(2) 



the coupling of the qubits to the environment. The 
strength of the coupling of qubit i is controlled by the 
parameter A^, while A = max{|Ai|}. 

We consider two different settings. In the first one 
acts on the space Hi H c , i-e. all qubits interact 
with a common environment. We shall name this the 
joint environment configuration. In the second one each 
qubit interacts with a separate environment. Thus the 
environment is split into n parts, H c = ®<=i %,t and 
VW, in eq. acts only on Hi ® H e ,i- This we call the 
separate environment configuration. The first case would 
be typical for a quantum computer, where all qubits are 
close two each other, while the second would apply to a 
non-local quantum network. 

We choose the initial state to be the product of a pure 
state of the central system (the QM) and a pure state of 
the environment. It may be written as 

|V>o) = \lpqm)\lpe), \lpqm) € H qm , \lp e ) G H c . (3) 

Thus, the reduced initial state in the quantum memory 
go = tr c \ipo)(ip \ = \ip qm }(ipqm\ is pure. In the separate 
environment configuration we furthermore assume that 
IV-'c) = ®i IV'e.i) with \ip e ,i) G H e ,i, which corresponds to 
the absence of quantum correlations among the different 
environments. 

For reasons of analytic simplicity we chose the purity 
P(t) as measure for decoherence. It is defined as 



P(t) = Tr pU(t), Pqm (t) = Tr e e- lllt \ii ){<ij; \e 



(4) 



To calculate purity or any other measure of decoherence 
we can replace the forward time evolution exp(-iHt) by 
the echo operator M(t) = exp(tHot) exp(—%Ht) where 



Hq determines the evolution without coupling (A = 0). 
In linear response approximation this operator reads as 



M(t) = l-i\I(t) - X 2 J(t), 



(5) 



where I(t) = f c \&tV t and J(t) = £ dr /J" dT'V T V T >. 
AVf = exp(iHot) XV exp(— iH^t) is the coupling at time t 
in the interaction picture. As discussed in [6| it is conve- 
nient to introduce the form p[p\ ® P2] = Tr(Tr c p\ Tr c p^). 
Purity is then given by 

P(t) = l-2A 2 f dr f dr'Rev4(T,T') + 0(A 3 ) (6) 
Jo Jo 

with 

A(r, t') = p[V T V T >Qo ® Qo] - p[Vr'goVr ® Qo] 

+ P[V T Q <g> Vr'Qo] ~ p[V T >Q0 ® Q V T ] . (7) 

Note that the linear terms vanish after tracing out 
the environment. Considering the form of the cou- 
pling in Eq. @ we can decompose X 2 A(t 1 t') — 
XiXjA^Xr, t') with 

A^\t,t') = p[V^V}P Qo ® Q ] - P [V^g V T ^ ® Qi ] 
+ P [V^g <8> V^g Q ] - p[V^ Qo ® QoVP] . (8) 

V t = exp(iH t) yW exp(— iH Q t) in analogy with V t . 
Equation ([6]) is then given as a double sum in the indices 
of the qubits. In the diagonal approximation (i = j), 
P(t) can be expressed in terms of purities P^{t) asso- 
ciated with individual qubits 



■it 

Pit) = i-E 



(9) 



P ( W («) = 1-2A? J dr J dr'A^(r,r') + 0(Xl). 

Each contribution corresponds to the purity decay of the 
central system in a spectator configuration with qubit 
i interacting with the environment. Eq. ^ is the cen- 
tral result of this paper. The diagonal approximation is 
justified in two situations: First if the couplings Vw of 
the individual qubits are independent from the outset, 
as would be typical for the separate environment config- 
uration or for the random matrix model of decoherence 
[6|]. Second, if the couplings in the interaction picture 
become independent due to mixing properties of the en- 
vironment, as would be typical for a "quantum chaotic" 
environment. 

To illustrate the case of independent couplings men- 
tioned above, random matrix theory provides a handy 
example. Such models were discussed in [^, 0, [l6| and 
describe the couplings VW by independent random ma- 
trices, chosen from the classical ensembles (l7j . 
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FIG. 1: (Color online) Purity decay for a GHZ state is shown 
for the environment described by eq. (113[l . The couplings 
are: all qubits to one spin (red triangles), nearby spins (blue 
stars), and maximally separated spins (green squares). The 
theoretical result (thick line) is calculated in terms of the 
Pup)(t) (thin lines). The main figure corresponds to the case 
of a mixing environment and the inset for an integrable one. 



In 0, Q the purity decay was computed in linear re- 
sponse approximation for two qubits, one of them being 
the spectator. For the sake of simplicity we choose the 
joint environment configuration, no internal dynamics for 
the qubits, and let H c and V^ 1 -* be typical members of 
the Gaussian unitary ensemble. Purity is then given by 



P ( W(t) = l-Af(2- Pl )/(t); 



(10) 



f(t) = t max{t,T H } + — min{t,T H } 3 . (11) 
6th 

Here, tr is the Heisenberg time of the environment, and 
Pi is the initial purity of the first qubit alone, which 
measures its entanglement with the rest of the QM. As 
we required independence of the couplings, we can now 
insert eq. (TTUj) in eq. © to obtain the simple expression 



P(t) 



(12) 



where pi is the initial purity of qubit i. In the 
presence of internal dynamics the spectator result is 
also known Q and can be inserted. As an example, 
we apply the above equation to an initial GHZ state 
(Q0---0) + |l---l))/\/^). Then all Pl = 1/2 and we 
obtain P(t) = 1 - (3/2)/(t) £V Af . For a W state 
(| 10 . . . 0) + 101 . . . 0) + • • ■ + 1 00 . . . l))/Vri, purity for each 
qubit is pi = (n 2 — 2n + 2)/n 2 m 1, in the large n limit, 
and purity decays as P(t) = 1 — f(t)Yli ■ For n = 2 
we retrieve the results published in [6| . 

The main assumption in Eq. ([9]) is the fast decay of 
correlations for couplings of different qubits to the bath. 
For the random matrix model discussed above this is 
trivially fulfilled. Yet, integrable environments are com- 
monly used [l8| and one may wonder whether Eq. ^ 



can hold in such a context. We shall therefore study a 
dynamical model where a few qubits are coupled to an 
environment represented by a kicked Ising spin chain us- 
ing identical coupling operators for all qubits. In this 
model, the variation of the angle of the external kicking 
field allows the transition from a "quantum chaotic" to 
an integrable Hamiltonian for the environment (l3j . The 
Hamiltonian of the chain is given by 



L-1 



L-1 



H (c) = ]T o^a^ + S(t) J2 b ■ u 



(e,i) 



(13) 



i=0 



where S(t) is a train of Dirac delta functions at in- 
teger times, L is the number of spins in the environ- 
ment, er( c ' 1 ) = (cri ' 4 -*, a y , cri ) are the Pauli matrices 
of spin i, and b the dimensionless magnetic field with 
which the chain is kicked (H = 1). We close the ring re- 



quiring a^'^ = cr (c ' 0) . The Hamiltonian of the QM is 

H (qra) = ^ fe . ^qm,^ where ff (qm,,) j s defined sim _ 

ilarly as for the environment. The coupling is given by 
= XJ2i ci qm ' i Vi c '" ,l \ where the ji's define the posi- 
tions where the qubits are coupled to the spin chain, e.g. 
if all ji are equal, all qubits of our QM are coupled to the 
same position. 

To implement a chaotic environment, we use b = 
(0.9,0.9,0). For these values the quantum-chaos prop- 
erties have been tested extensively 14|. We use a ring 
consisting of L — 12 spins for the environment and 4 
additional spins for the qubits of the QM. The coupling 
strength has been fixed to A = 0.005. In Fig. [TJ we 
study purity decay when all four qubits are coupled to the 
same spin, when they are coupled to neighboring spins, 
and when they are coupled to maximally separated spins, 
j — (0, 3, 6, 9). The initial state is the product of a GHZ 
state in the QM and a random pure state in the envi- 
ronment. We compare the results with the prediction of 
Eq. ([9]) obtained from simulations of the corresponding 
spectator configurations (thin solid line). All spectator 
configurations yield almost the same result for P(t), so 
that we can see only one line. The figure demonstrates 
the validity of eq. © for well separated and hence in- 
dependent couplings. For smaller separations we obtain 
faster decay. 

Next we perform similar calculations for integrable en- 
vironments. A typical case is shown in the inset of fig. [T] 
for a two qubit QM and b = (0,1.53,0). Both qubits 
are coupled to the same, neighboring or opposite spins in 
the ring. Purity decay is faster than in the mixing case 
as can be expected from general considerations in [l9j |. 
Nevertheless the sum rule is again well fulfilled except if 
we couple both qubits to the same spin. This leads us to 
check the behavior of the correlation function 

P[V^V}P Q0 ® go] = (MVVvPlfh). (14) 

This is the simplest and usually largest term in AW, 
eq. ([5]). Figure [5] shows this quantity for mixing and 



4 




20 40 20 40 20 40 

T 



FIG. 2: (Color online) The color-coded absolute value of the 

real part of {4>o\Vt^ V^P \i/jo) given in eq. (I14|l as a function 
of r and r' for the integrable [(a) (b) (c)] and the chaotic 
[(d) (e) (f)] spin chain environment (L = 8). (a),(d) show the 
autocorrelation function (i = j); (b), (e) and (c) (f) the cross 
correlation function for couplings to the same spin or opposite 
spins in the ring respectively. 

integrable environments in the first and second row re- 
spectively. The first column shows the autocorrelation 
function (i = j). The second and third column give the 
cross correlation function (i ^ j) when the qubits are 
coupled to the same or opposite spins respectively. We 
see that in the latter case correlations in both the inte- 
grable and the chaotic case are small, thus showing that 
also in integrable situations our condition can be met. 
Non-vanishing correlations, as in fig. EJb) and fig. [2]Je) 
lead to deviations from the sum rule. The pronounced 
structure for the integrable environment, fig. [DJe), may 
be associated with the oscillations of purity decay (inset 
of fig. [T|). 

We have calculated generic decoherence of an n-qubit 
quantum memory as represented by purity decay. For 
large purities it is globally given in linear response ap- 
proximation by a sum rule in terms of the purities of 
the individual qubits entangled with the remaining reg- 
ister as spectator. This sum rule depends crucially on 
the absence or rapid decay of correlations between the 
couplings of different qubits in the interaction picture. 
We prove, that this is fulfilled by sufficiently chaotic en- 
vironments or couplings, but using a spin chain model as 
an environment we find, that even if the latter is inte- 
grable correlations can be absent and the sum rule holds. 
While exceptions, no doubt, can be constructed we have 
a very general tool to reduce the decoherence of a QM of 
n qubits to the problem of a single qubit entangled with 
the rest of the QM. Furthermore the result is equally 
valid if we perform local operations on the qubits. The 
sum rule also applies if the register is split into arbitrary 



sets of qubits. In fact one and two qubit gates are known 
to be universal for quantum computation [l|, so we lay 
the foundation for the computation of decoherence dur- 
ing the execution of a general algorithm. This extension 
only requires the knowledge of the decoherence of a pair 
suffering the gate operation while entangled with the rest 
of the register. As each gate is different, we will have to 
take a step by step approach, which at each step will 
only involve one and two qubit decoherence. The linear 
response approximation will be sufficient due to the high 
fidelity requirements of quantum computation. 
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